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Abstract 

Borel summable semiclassical expansions in ID quantum mechanics are considered. 
These are the Borel summable expansions of fundamental solutions and of quantities 
constructed with their help. An expansion, called topological, is constructed for the cor- 
responding Borel functions. Its main property is to order the singularity structure of the 
Borel plane in a hierarchical way by an increasing complexity of this structure starting 
from the analytic one. This allows us to study the Borel plane singularity structure in a 
systematic way. Examples of such structures are considered for linear, harmonic and an- 
harmonic potentials. Together with the best approximation provided by the semiclassical 
series the exponentially small contribution completing the approximation are considered. 
A natural method of constructing such an exponential asymptotics relied on the Borel 
plane singularity structures provided by the topological expansion is developed. The 
method is used to form the semiclassical series including exponential contributions for 
the energy levels of the anharmonic oscillator. 
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1 Introduction 



In this paper we continue our investigations to represent basic quantities of the quantum 
mechanics in the form of Balian - Bloch representation i.e. in the form of the Laplace-Borel 
transforms in which the cojugate variables are an action and the Planck constant 0. The 
key results, which the paper is relied on, has been published earlier Q. The present work 
develops these key ideas and using the explicit form of the fundamental solutions [|ll], 12] to 
the ID Schrodinger equation expresses the Balian-Bloch representation in the form of what 
we call a topological expansion. We describe also the way of using the representation to 
construct extended JWKB approximations in the form of so called exponential asymptotics 
(sometimes called also the hyperasymptotics, see [23, 24, 25, 26, 27] and the references cited 
there) and we consider some particular applications of the Balian-Bloch method as well. 

For simplicity the potentials considered in this paper are assumed to be polynomial but 
its main results are valid for more general meromorphic potentials as well. 

Being Borel summable the fundamental solutions define their corresponding Borel func- 
tions with the help of which they can be represented in the form of the Borel transformation 
from the Borel plane of the action variable to the complex plane of the ft _1 -variable. For the 
polynomial potentials these Borel functions are in fact all the same despite the fact that they 
are defined by different fundamental solutions ||. This means of course that the fundamen- 
tal solutions themselves are in close relations to each other being in fact a mutual analytical 
continuation of each other in the ftT 1 plane [||, [|. 

Therefore, to get any of the fundamental solutions it is necessary only to choose properly 
an integration path in the Borel plane. However, to do it a detailed knowledge of singularity 
distribution of the Borel function in the Borel plane is necessary. It is the aim of this paper 
to provide us with an effective tool for studying these singularities. Namely, we develope an 
expansion for the Borel function called toplogical in which an expansion parameter is the 
complexity of the Borel plane corresponding to successive terms of the expansion. 

With the help of the fundamental solutions we can solve most of the ID quantum me- 
chanical problems so that the corresponding quantities involved in the problems considered 
depend on different pieces of the fundamental solutions used. These quantities themselves 
can have then semiclassical representations which can be Borel summable and can serve as a 
source of their semiclassical approximations as well. It is clear that the corresponding Borel 
plane singularities of these quantities are defined then by the pieces of the fundamental solu- 
tions constructing them. Therefore the topological expansion method can be applied also to 
determine the approximate singularity structure for these quantities as well. 

The semiclassical expansions used as a source of approximations are considered as insuf- 
ficient providing us with unavoidable nonvanishing errors. It is well known that the reasons 
for these errors are immanent due to the divergence of the semiclassical series so that the 
latter as asymptotic neglect the exponentially small contributions. Nevertheless, since the 
series are Borel summable they have to contain the full information about such exponential 
contributions. A common goal of many approaches was just to recover these contributions 
leading to a formulation of so called resurgent theory |25|, 26, 27]. 

Let us note, however, that the exponentially small contributions is of its own importance 
since in many cases of quantities considered these contributions are dominant. Among the 
latter cases the most well known one is the difference between the energy levels of different 
parities in the symmetric double well p0| . But these are also the cases of transition probabil- 



ities in the tunnelling phenomena 3C] or their adiabatic limits in the time dependent problem 
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of transitions between two (or more) energy levels (see [[32|, |33| and references cited there) or 



the exponential decaying of resonances in the week electric field (see [34, 35 1 and references 
cited there). 

In an approach of our paper we make full use of the Borel summability of the quantities 
considered as well as of the corresponding topological expansions to construct the relevant 
exponential asymptotics. 

However, as a necessary step of our fomulation is the knowledge of the Borel plane singu- 
larity structure of any considered quantity. It is just the topological expansion which allows 
us to built this knowledge step by step. 

The toplogical expansion is constructed directly from the Froman and Froman represen- 
tation of the fundamental solutions which themselves are given in the forms of functional 
series |2[ |3[ Q. Therefore, we shall start in the next section with the detailed presentation of 
the series. 

In Section 3 the topological series representation for the Borel functions is introduced and 
its convergence is proved. This representation provides us with an algorithm for approximate 
calculations of Borel functions being alternative to the ones relied on the Pade approximants 
P, B pi, [j~4| ], continued fractions or conformal transformations [14|. 



In Section 4 singularity structures of the topological series expansion are analyzed and 
their hierarchic form (which gives rise to the name of the series) is established. We consider 
there as the simplest examples the 'first sheet' singularity structures of the linear and har- 
monic potentials. In particular we describe completely the singularity structure of the Borel 
plane of the harmonic oscillator Joos function found first by Voros by a different method [^] . 

The results obtained in Sections 3 and 4 are applied in the next section where the solution 
of the so called connection problem within the framework of the Balian-Bloch representation 
is discussed. 



In Section 6 we discuss the problem of the exponential asymptotics ]23| , 24, |2^, E7|]. 
We show that this problem has a natural solution in the framework of the Balian-Bloch 
representation and gets a natural support from the topological expansion method. 

In Section 7 the energy levels of the single-well anharmonic potential are considered in 
order to show how to use the topological expansion to construct their extended exponential 
asymptotics. 

Finally, in Section 8 we summarise our results. 

2 Fundamental solutions 

Let us remind basic notions of our considerations (see for details). 
The fundamental solutions satisfy the Schrodinger equation: 

V"(x, A, E) - X 2 q{x, E)V(x, A, E) = (2.1) 

where: q(x,E) = V(x) — E, A = \j2mftT 1 . Both A and E can take on complex values. V(x) 



is assumed to be a polynomial of any degree n > 1. A Stokes line pattern (see [11, [T^] for 
necessary definitions) relevant for our considerations is shown in Fig. 1. (A total number p 
of sectors equals to n + 2 in this case.) 

The following fundamental solution ^f(x, A, E) to ( ^l| ) can be attached to the sector S%: 

n (x, A, E) = q-Hx, E)e° X C qh iy ' E)dv xf(x, A, E) 
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o" / q 2 {y,E)dy 

J Xo 



< 0, x & Si, a = ±l 



(2.2) 



q(x ,E)=0 

with the "amplitude factor" Xi( x i given by the following functional series: 



Xl(x, X,E) = 1 + ( 2^ J J dyi... \ dy n u(y x ) . . . u(y n ) 

. \ _ e 2AC( yii x)l . U _ e 2\aV2,yi) 



(2.3) 



1 -e 



2A5(j/„,j/„_i) 



where: 



q"(y) 5g%) 



'*(!/) (« *(!/))' 



(2.4) 



£(zo,aO = -o- / qi(y,E)dy 

and where an obvious dependence of u, q, £, etc. on E has been dropped. We shall also put 
a = — 1 in ( |2.2[) -( ^4|) assuming that in (|2.2|) the corresponding inequality is satisfied in this 
case. 




Rex 



Fig. 1. The Stokes graph for a general polynomial potential 
In the sector S\ the following semiclassical expansion for Xi( x i ^) takes place: 

«l,n(aj) 



Xl (x,A)~ X f(^A) = l + ^ 



n> 



^ 2A)» 
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«i,n(a?) 



dUq * (6) |5 * (&> 



00 fc 



(2.5) 
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fe = l,2, 



As it has been shown in if x stays in S\ of Fig. 1 then we can define for 9?s < the 
following Laplace transformation of the amplitude factor xi(x,A): 



(2.6) 



with the integration contour C shown in Fig. 2. (The factor 2 in the exponential function in 
( |2.5| ) is introduced for convenience). By the form ( |2.6| ) Xi(x, s) is defined holomorphically in 
the half-plane "Sis < $1£(xq,x) and since $l£(xo,x) is positive x\(x,s) appears to be, in fact, 
the Borel transform of A). The contour C in (^1]) can be chosen as a circle with its 
radius A to be large enough to substitute xi(a;, A) by its semiclassical series fl2.5|) . Then for 
\s\ < |^(xo,x)| the LHS of (|2.6|) can be integrated to give the following Borel series: 



(- s y 



n>l 



(2.7) 



convergent in the circle |s| < |5ft£(a;o, x)\. The point sq(x) = £(xq,x) is a singularity for 
Xi(x,s) closest to the origin. 

Im X 
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Fig. 2. The cut A-plane corresponding to the global solution ^f(x, A) 
The transformations (|2.5| ) can be inverted to give: 



Xi(x,A) = 2A le ZAs x 1 (x,s)ds 
Jc 



where the contour C starts at the infinity 3?(As) = — oo and ends at s = 0. Since C can be 
freely deformed in the half plane 5Rs < the formula (2.8) define Xi( x i A) in the whole sheet 
shown in Fig. 2 excluding the points of the negative real half-axis. 
The following is worth to note. 

The formula (2.8) is certainly valid for 5RA > when the contour C stays in the half 
plane Ifts < 0. It can be continued, however, to other domains of the Riemann A-surface 
corresponding to Xi( x > A) if accompanied with suitable changes of the variable x. Thus, for 
example, when continuing x to the sector S2 and deforming the contour C in Fig. 2 into C\ 
the formula ( |2.5D will then define Xi(x, s ) m the half plane > 0. On the other hand the 
inverse formula (2.8) defines then Xi( x i A) in the half plane 5RA < with the contour C in the 
formula deformed (anticlockwise) from its position in the left half-plane to its new position 
in the right half of the s-plane. The function Xl{ x i^) fulfils then for A > the condition: 
Xi(x,— A) = X2{x,\). Possible singularities of xi(x, s) existing in the corresponding half- 
plane > when x stays in the sector Si move to the half-planes 5Rs < when x moves to 
the sector £2 (see also Section 5 for a relevant discussion). 



3 Topological expansion of Borel function xi(x,s) 

As it follows from the definition of xi(x,s) if we want to learn something about it we 



have to analyze Xi( x > A) as given by (2.3). We shall show below that if x stays in S\ (see Fig. 



2) then it is possible to represent each term of the series in (2.3) in the form of the Laplace 
transformation i.e. we shall show that Xi( x i s ) can also be represented by some convergent 
functional series. The series however can still be continued to almost the whole x-plane when 
the latter is deprived of some vicinities of their turning points. 

To this end let us consider the n th term of the series in fl2.3|) and particularly its re-fold 



integral. Introducing to it £ = = £(xq,x) as a new integration variable we get: 



Y n (t,\)= J #1... J de n ^ 1 )...^ n )(i- e 2A ^))...(i- e 2A ««"- 1 -^) (3.1) 

7l(f) 7l(£n-l) 

where u(£(x)) = uj(x)q~^ (x) and the path 7i(£) starts from = +00 and ends at the point 
£ of the £-plane. 

Let us point to some basic properties of the transformation £ = £(#). First it maps the 
two sheeted Riemann surface which the x-plane actually is into another (in general infinitely 
sheeted) Riemann surface on which each sector of Fig. 2 is represented by the right (left) half 
planes. In particular sector 1 in Fig . 2 is mapped into the right half of this cut £-plane whilst 
sectors 2, 3 and p into the left ones (see Fig. 3 where the sectors 2 and p lie below the sheet 
shown). Zeros of q(x) which are singular points for uj{x)q~^{x) are also suitably transformed 
into the corresponding root branch points (of the third degree) of cD(£(z)) on the £- Riemann 
surface (see Fig. 3 and Appendix 2). On this surface a)(£(x)) becomes additionally infinitely 
periodic with its complex periods acting however between different sheets of the surface. As 
a result of this an image of each root of q(x) proliferates infinitely on the ^-Riemann surface 
with every such a copy giving rise to still new branch point and sheet. The only exception 
of the latter rule is the linear potential case the ^-Riemann surface of which is three sheeted 
with a single root branch point of the third degree. 
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Fig. 3. The £-plane singularities corresponding to <!-[ (£, s) (case q = 0) 



Opening the brackets in (JO) we get: 



^«(£,A) = £ ^;n...r 28 (e,A)(-l) 

< T\ < . . . < n 
0<2q<n 



ri-r2+r 3 -...+r 2q -i-r 2q 



where 



with 



and 



71 (0 7l(5n-l) 



~1 I 



o(0= y dTjwfa) 

71 (0 



*W..r 29 (£,A) = | #1 ■ ■ ■ / ^(6) • • • ^n)^" 1 -^+^3-"+^,-i " 



71 (0 7l(Sn-l) 



9 = 1,2,3,... 



Note, that all the integrals in ( |3.5|) a re absolutely convergent. Therefore, it should be 
now obvious that to each integral in fl3.5|) the following Laplace transformation form can be 
given: 



*W..r 2 ,(£,A) = J dse 2Xs Y n . ri ... r2q (C,s) (3.6) 
c 

where the integration contour C starts at $ts = — oo and ends at s = and the Laplace 
transform Y n . ri __ T2q (^,s) is to be determined. We do it in Appendix 1. An important ob- 
servation done there is that it is possible to rearrange the order of terms in the series (|2.3| ) 
in such a way to sum it in an accordance with the increasing q rather than n - the number 
of the integrations in ([^). (All these are still possible since the series ( [2.3| ) is absolutely 
convergent). As a result of such reordering Xl(d A) can be represented as the following sum: 

Xi(C,A) = 2A^xi' ?) (C,A) (3.7) 

q>0 

£ G Si, | arg A| < tt 

where 

X?(£,A) = Jdse 2Xs <b[ q \n,s) (3.8) 



with ¥ q \£,s), q > given by formulae (|A1.12|) of Appendix 1 and with the contour C shown 



in Fig. 5. Of course, since the series ( |3.7D is absolutely and uniformly convergent we have 
also: 



Xi(£,A) = 2A / dse 2Xs ^,s) (3.9) 



c 



with <E>i(£,s) given by ( |Al.ll|) so that the corresponding Laplace transform Xi(£> s ) defined 



by ( |2.q ) can be identified as: 

Xi^,s) = ^,s) (3.10) 

The expansions ( |3.7| ) and ( |A1.10|) shall be called further topological expansions for the 
following two reasons: 



1. the higher term of the series in ( A1.12Q , the more complicated is its Riemann surface; 

2. the Riemann surface R q corresponding to the term s) in (ALII) can be reduced 



to some R q ' with q' < q when deprived of some singular points of (£, s) i.e. a set 

S q of all singularities of $^(£, s) includes a set S q > corresponding to ^^(S,, s) (see the 
next section). 

3.2. Analytic properties ofxi(£,s) 

The analytic properties of Xi(£) s ) have been established in Section 3 of Appendix 1. As it 
follows from Appendix A1.3 the Laplace transform xi(£) s ) is holomorphic in some vicinity of 
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the point s = for £ E R(tf') i.e. it is the Borel function Q2.7|) corresponding to Xi(£> A). F° r 
5R£ > 0, however, xi(£, s) is holomorphic in the half plane < 0. Therefore, the asymptotic 
series constructed for Xl(£> A) when A —* oo is Borel summable to the function itself - a result 
which is in a full accordance with the corresponding one obtained in P] and mentioned in 
Section 2. 



4 Singularity structure of s) 

Because of (|3.1p) this is the singularity structure of <3?i(£,s) and the latter structure is 
determined by the corresponding singularity structures of (£, s) due to (|A1.10|) . These 
structures can be investigated by the analytic continuation procedure of the formulae ( |A1.11| ) 
- ( A1.12j ) with respect to s and £ and are, on its own, determined completely by the corre- 



sponding singularity structures of and the integrations present in (Al.ll) and ( |A1.12| ) 



(see Appendix 1). These integrations can give rise to singularities due to the following two 
mechanisms [Q|: 

1. moving singularity of the integrand approaches a fixed limit of the integration or, in- 
versely, a moving limit of an integration approaches a fixed singularity of the integrand 
(so called end point (EP-) singularities). 

2. moving singularity of the integrand approaches some another singularity pinching un- 
avoidably in that way the integration contour (so called pinch (P-) singularities). 



In the convolution integrals of the formula ( Al.ll ) only the functions <!>(£) and can 
give rise to both the (EP- and P-) singularity mechanisms since a dependence of the integrals 
on the remaining partners of the convolutions are holomorphic. 



From the defining formulae ( [A1.12 ) and from the £-Riemann surface structure on which 



cj(£) and are defined (this structure was sketched in the previous section) it follows also 
that even for the simplest cases of first few (£, s)'s their global (£, s)-Riemann surface 
structures are too complicated to be fully handled and only some crude descriptions of them 
are possible limited to a few first sheets and a few singularities. 

However, in making the corresponding analysis by limiting ourselves to first few q's we 
are free in deforming the integration contours in ( A1.12| ) i.e. the limitation of 7i(£) to the 



canonical choices is no longer valid. This observation is very important and prooves that 
the Borel function Xi(£, s) constructing initially for the fundamental solution of the sector 



Si is universal i.e. each Borel summable solution to the Schrodinger equation (2.1) can be 
obtained by the Borel transformation of Xi(£, s) with a properly chosen integration path in 
the Borel plane. A discussion of the latter property of the Borel summable solutions and 
some of its consequences is postponed however to another paper ||. 

Having in mind the incredible (in general) complexity of the (£, s)-Riemann surface struc- 
ture of <E>i(£, s) we shall describe first a general procedure of getting this structure for first few 
i>i (£, s)'s taking into account also a few singularities of u)(£) and fi(£) and next we try to 
give as full as possible a description of such structures for the linear and harmonic potentials. 

q = 
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It is seen from (|A1.12j ) that ^(^s) is an entire function of s for any ^ not coinciding 
with singularities of u)(£). Its singularities in the ^-variable coincide therefore with those of 
0(£) and consequently with those of cD(£) as the EP- singularities shown in Fig. 3. 

q = i 

This is the Riemann surface structure of s) as defined by ( |A1,12| ) for q=l. 

~m, , r , w Ji([8( 8 -77)nK-»7)-4( a -77)n(0]5) 

$ (Z,s)= dr](j^- V )(2s-2ri)^ ^ 4.1 

Jc(s) [8(s- V )Q(C-v)-^s- V )n(0}^ 

It follows from ( |4.1| ) that singularities of the subintegral function are essential singularities 
coniciding with the branch points of fi(£) and f2(£ — if). Because of the single rj- integration 
in (4.1) only the EP-mechanism can generate singularities in the 's-plane' since all the rj- 



singularities are the moving ones (depending linearly on £) so that the positions of all the 
(essential) singularities of $^ (£, s) coincide again with those of 0(£ — s) and f2(£). Therefore, 
these positions on the £, s-Riemann surface are the following: 



etc. 



(4.2) 



The nature of all these singularities is not altered by the integrations i.e. all they are 
branch points. Therefore, the resulting pattern of cuts on the corresponding Riemann surface 
which follows from Fig. 3 is sketched in the figuers 4-5. 

q = 2 

This is the Riemann surface structure of &\ (£, s) as defined by (Al.ll) for q = 2. From 



(A1.14) we have: 



(£,*) = J d V J d^ 1 -rj)u^ 1 )(2s-2 V ) 2 ^^- (4.3) 

C(s) 7(5) 

z = 4(s- rj) (fi(£) - 20(6) + 20(6 - V)) 

Note that the ^-integration in ( [4.3[ ) runs across a sheet of the £-R<i emann surface shown 
in Fig. 4 (where the s variable is to be substituted by the rj one). However, contrary to the 
close correspondence between the distributions of sectors and turning points on the Stokes 
graph of Fig. 1 and of sheets and the corresponding cuts on Fig. 3 such a correspondence is 
lost in the case of Fig. 4 i.e. we are left only with some properly arranged system of branch 
points and cuts. 

Since the half of the cuts in Fig. 4 are moving then except of the EP-singularities the 
P-singularities are also generated by both the £- and reintegrations in ( |4.3| ). 
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s+C 2 Jc 




Re^ 



Fig. 4. The £-plane singularities corresponding to $1 (£, s) (case g = 1 



Consider first results of the ^-integration in ( |4.3| ). 

The EP-singularities which follow from this integration coincide (with the corresponding 
substitution s by 77) with those in the figuers 4 and 5 are given again by (|4.2|). 

* Ims 



/ \ 




Res 



Fig. 5. 



The s-plane singularities corresponding to l>^(£, s) (case q = 0) 



A generation of P-singularities can be performed by moving singularities depending on 77 
(see Fig. 4). For example, moving clockwise the singularity 77 + Ci around the end point £ 
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of 7i(£) and next pinching 7i(£) against £1 we generate a singularity of ( ]4.3| ) at r\ = in the 
77-Riemann surface. It is placed however on another sheet of the surface since to achieve it 
we had to go around the branch point singularity £ — Ci> shown in Fig. 5, in the clockwise 
direction. 

To obtain all other 77-plane singularities generated by the ^-integration in ( |4.3| ) we proceed 
in the same way as described above. All these singularities lie on sheats which can be reached 
by going around the two branch points (in any direction - clockwise or anticlockwise) shown 
in Fig. 5. Therefore all these singularities are shared by the actual positions of the branch 
points cuts of Fig. 5. They can become visible by cutting the 77-plane in a different way or 
moving appropriately both the branch points to the left. 

Im % 




Fig. 6. The £-plane singularities corresponding to subintegral function in (4.3) 



Choosing for example the last possibility and moving £ toward Sector 3 we shall arrive 
at the situation shown in Fig. 6. If £ and 77 are moved so that 3i£ < = %1(t) + Ci) then a 
further motion of 77 + Cl upwards to the point (2 pinches the path 7i(£) producing in that way 
a singularity at 77 = £21 — C2 — Ci- It lies to the right from the cut at £ — £1 in the '77-plane' 
and is screened therefore by the cut just mentioned when 3i£ > KC2 (see Fig. 7). 

By the identical analyses applied to each pair 77 — Q, Q of the singularities lying on the 
sheet in Fig. 4 the singularities at s = Cij or at s = Cji = ~dj can De produced being screened 
by cuts at s = £ — Cj or at s = £ — Ci, correspondingly. All the singularities produced in this 
way are branch points. 

According to ( |4.3j ) the second, final integration is performed over the 77-plane providing 
<3?i (£, s) with all its £- and s-plane singularities. This integration transforms all the 77- 
singularities obtained by the first (£-) integration into the corresponding s-ones (by the EP- 
mechanism) and provides us with additional ^-singularities by the pinch mechanism. Pinching 
for example the singularity £ — C2 against (21 we obtain the ^-singularity at £ = (2 + C21 lying 
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on a sheet of the £-Riemann surface originated by the branch point at C2 on Fig. 4. This 
branch point is screened, of course, by the cut at £ = s + C2 when IRs > K(y (see Fig. 8). 
Therefore, figures 7-8 show the complete singularity structure of Q\ (£, s) when continued in 
£ in the way shown in Fig. 6. 




Im s 



Res 



~ (2) 

Fig. 7. The s-plane singularities corresponding to &\ (£, s) 



^.i. T/ie analytic structure of the Borel function for the linear potential 

We can put for this case q(x, E) = x and £ = x 3//2 with the corresponding Stokes graph 
shown in Fig. 9. and we shall consider as the Borel function defined by the funda- 

mental solution ^?i(x, A). 

Im % 



s+£ ? 




Fig. 8. The £-plane singularities corresponding to (£, s 
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At the first glance the corresponding analysis seems to be simple because of the simplicity 

of the relevant functions = — jQ~gz and ^(C) = | as a result of which the three sheeted 

Riemann surface branching at £ = (the surface being the image of the two sheeted x-plane 

by the transformation £ = x 3 / 2 ) decouples into three independent sheets. The unity of the 

surface is recovered however by the solution ^i(x,\) which being holomorphic at x = 

branches at this point as £ 2 / 3 when considered as a function of £. However, since we are 

interested in the properties of the Borel function <I>i(£,s) determined rather by xi(£> A)/2A 

it is the latter the (£, A)-dependence of which is most important. 

* Imx 




Rex 



Fig. 9. The Stokes graph for the linear potential 

The latter dependence can be established to some extent noticing that continuing analit- 
ically the solution V&(£, A) = £ ^e~ A £xi(£i A) in the A-plane (whilst £ is fixed) by rotating A 
by the angle ±6tt we come back with the beginning of the integration path 7(£) in A) 
to the infinity of the first sector. Of course, this path is by the above continuation deformed 
from the initial canonical one into the one surrounding the point £ = twice (in the direction 
suitable to the sign) to end eventually at the point £. This is because continuing Xi(£;A) 
in A in the above way we have to shift the infinite end of the path to the neighbour sectors 
each time when A changes by ±7r (this operation keeps the factor e~ x ^ always vanishing in 
the infinities of the passed sectors) @ . 

However, the above A-continuation of Xi(£) A) is equivalent to its continuation to the same 
point £ along the deformed path starting from its initial canonical form. Since by this 
latter continuation the argument of £ changes also by ±6ir then the factor £~6 of ^(^A) 
aquires minus by this continuation so does the factor Xi(C)A) since by this continuation 
A) can not change because it branches at £ = as £ 2 / 3 . It follows therefore that Xi(£i A) 
branches at £ = as £ 1//6 . 

From the latter observation it follows further directly that the Borel function <I>i(£,s)(= 
Xi(£,s)) branches at the infinity point of its s-plane also as s 1 / 6 . This can be seen noticing 
that to recover the factor xi(£; A) by the Borel transformation of < I ) i(£, s) we have to change 
successively the integration path in the transformation from the negative real halfaxis to the 
positive one (and vice versa) according to which the sector the infinite end of the deformed 
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path 7(£) is actually in. These Borel transformation paths are again the deformations of each 
other obtained by moving the infinite end of them along the circle of infinite radius i.e. all 
the singularities of <&i(£,s) are avoided by these deformations. Since after six such changes 
the Borel transformation of s) has to change its sign in comparison with its initial value 
so $i(£, s) itself has to do it. 

Therefore, we conclude that for fixed £ the s-Riemann surface of <3?i(£,s) is built of six 
sheets. 

The above situation is however not so simple when the formulae ( A1.12[ )-( |A1.14 ) defining 
s) are considered. The Bessel functions in these formulae convert the simple pole of 
0(£) at £ = into a corresponding root (of the forth order) branch points accompanied by 
essential singularities (see Appendix 3.1). Also the successive £- and reintegrations in these 
formulae have to generate unavoidably the branch points at £ = 0, s = and £ = s of the 
logarithmitic type. This is of course because the representation of s) given by ( A1.12| )- 



(Al.14) is singular providing us with the correct positions of singularities but not necesserilly 
with their nature. The above example of the linear oscillator shows that the proper behaviour 
of <3?i(£,s) close to its singularities is obtained only by the full resummation of these series. 
Nevertheless, in more complicated cases of potentials an information the series provide us are 
certainly very useful. Also in the case just considered. 

Namely, taking into account the reccurent relations ( |A1.14 ) we can establish inductively 
that 3>i(£, s) being defined on its six sheeted (£, s)-Riemann surface has on its first two sheets 
singularities shown in the figures 10a, b (see Appendix 3.1 for details). The point s = 
on the sheet of Fig. 10b is regular for <£i(£, s), according to general results of App. 1. 
According to this analysis the points £ — s = are the four order branch points of 3*1 (£, s) 
and simultaneously its essential singularities but we should have in mind that the last two 
properties can be incorrect. 

The same property concerns the points £ = and s = the latter being on the second 
and further sheets of Fig. 10b. All these points arrange themselves to build in the considered 
approximation of <l>i(£, s) infinitely sheeted Riemann surface. However, even for this simple 
case the topology of the surface except its first two sheets is too complicated to be fully 
described. 



Im \ 

7© 



Im s 



Re!; 



Res 



Fig. 10. The £- and s-plane singularities corresponding to $^(£, s 
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Nevertheless, one general conclusion valid at least for all the polynomial potentials can be 
drawn from the above consideration. Namely, if for a general polynomial potential we consider 
any pair of neighbour sectors joined by the analytic continuation in A when A — > e ±l7T A and 
we continue a fundamental solution defined in one of the sectors to the second along the 
canonical path then the corresponding 's-plane' singularity structure of the first sheet of the 
respective Borel function s) is exactly the same as for the 'simplest' case of the linear 
potential described above. 

4-2. An alternative non-standard Borel representation for the linear potential wave func- 
tion 

In the previous subsection we have made a disappointed note that even in such a simple 
case as the linear potential one the corresponding Borel function properties which follow from 
the topological expansion are quite complicated. We have however shown also that the actual 
structure of the linear potential Borel plane should be rather simple. Below, we want to show 
that indeed this complication is apparent and changing a little bit the definition of the Borel 
function one can simplified the latter enormously for the case considered. Namely, let us 
replace the definition ( |2.7| ) of the Borel function by the following one 



n>0 1 V 1 i " 2> 

which corresponds to the following representation of xf H (£, o~) by the Laplace transformation 



+ioo+Ao 



TTl J 



e - 2 A.Xl(i 1 A) dA 



— ioo+Ao 

< A Q < 1, 



(2A) 



(4.5) 



a < 



so that the invers Borel transformation is given by 



Xl (£,A) = (2A)§ / e 2 ^xf\i,a)da (4.6) 

J— oo 

Let us now make use of the fact that the fundamental solution x I / i(x, A) can be given the 
following integral representation (see p0|| , Mathematical appendix) 



*i(x,A) = -t=(2A)^ / e^-^dy (4.7) 

'TV Jc 



where we put x real and positive and the contour C is shown in Fig. 9. 

Changing in ( [4.7D the integration variable y into x~ l l i y and next putting 2a = x 3//4 y 
x _3//4 y 3 /3 + 2x 2 / 3 /3 we can bring the integral to the following form 
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*i(x,A) = x 



i 

"4e" 



„ 3 



(4.8) 



1 3 3 3 / ~ 2 3T\ 3 / 1 3 3 3 / ~ _ _ , 

-3(<7 — -x2 )xi + 3x4 y (7^(7 — _x 2 j I - I -3(cr - -X2JX4 - 3x*ya(a - -x?) \ da 
Hence for Xi lt (£.>°~) we get finally 



xf (M 



(-3(a - JO(iO* + 3(|0* V^O) 3 (4- 

-3(o- -k)(io»- 3(16* v^ef 



It follows from flO) that xf*(£, a) is denned on the two sheeted Riemann surface having 
the branch points a = and c = £ as its unique singularities. 

The non standard representation fl4.5|) -( [03l ) of the Borel function considered above shows 



that the complicated form (2.7) of the standard one depends on the representation itself and 



it can be simplified greatly by the proper choice of such a representation. 

4-3. The singularity structure of the Borel function for the harmonic oscillator 

Making, if necessary, a suitable rescaling we can put in this case q(x) = x 2 + 1 (assuming 
the energy to be negative). The corresponding Stokes graph is then shown in Fig. 11 and 
we choose as usually the sector 1 to provide us with the fundamental solution ^i(x, A) and 
its Borel function $i(£, s). Because of the last conclusion of Section we consider now a 
case of the Riemann surface structure corresponding to 3>i(£, s) when £(= Jf i y/y 2 + ldy) 
is continued to the sector 3 of Fig. 11 (along a canonical path). Then the first sheets of 
$[ 1} (e,s) and <±H 2) (£,s) are shown in the figures 12a,b and 13a, b respectively. 

Using again the formulae ( |A1.14| ) we can show inductively that the first sheets of $^ 2< ^ (£, s) 
and <li 29+1) (£,s) look as in the figures 14 and 15. All the detailed considerations establishing 
these can be found in Appendix 3.2. 

4-4 Borel plane structure of harmonic oscillator J 00s function 

When in the consideration of the previous subsection we shall push K£ to minus infinity 
(this corresponds to push x to the infinite point 003 of the sector 3 of Fig. 11) than we get the 
'Borel plane' singularity structure of the so called Joos function for the harmonic oscillator. 
The last name is given to the coefficient Xi^W = h m £^oo 3 ^) (HI so that the energy 
spectrum of the harmonic oscillator is given by xi—^A) = 0. Note that in the limit £ — ► 003 
all the functions $^ 29 (£, s) vanish so that the corresponding limiting functions &i^ 3 (s) 
contribute only to Xl-*3( s )- 
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Im x 



Y,(x) 



Rex 



Rex 
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Fig. 11. The Stokes graphs correspond- 
ing to the harmonic oscillator 



As it follows from the considerations of the previous subsection the singularity structure 
of the latter function is determind by the branch points distributed along the imaginary axes 
of the s-Riemann surface. This distribution can be discribed completely if instead of the 
Borel function Xi-+3{ s ) we sna l consider the one corresponding to logxi^3(A). To this end 
let us note that as it follows from Fig. 11a the normal sector of xi—^A) (i.e. the one where 
Xi-»3(A) is holomorphic and can be expanded into the semiclassical series (|2.5|)) is defined by 
| argA| < 7r. One can easily find also (by analytic continuation in A) that 



X2^ 4 (A) = Xi^s(Ae^), < | arg A| < tt 



(4.10) 



where a = argA/|argA| and XiL^A) are the canonical coefficients corresponding to the 



graphs of Fig. lib and 11c respectively. Despite (4.10) these two canonical coefficients obey 
the following two other relations 
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Xi-3(A)X2^ 4 (A) = 1 + e mX , < arg A < 7r (4.11) 

Xi->3(A)xt,4(A) = 1 + e"™\ -vr < arg A < 



in which the fact that 2 JV \/y 2 + Idy = Tri has been used. 

The relations ( |4. 11 ) follows as a result of an identity which the four fundamental solutions 



^k, k = 1, 4, corresponding to the Stokes graphs of Fig. 11 have to satisfy since only two 
of them can be linearly independent. 



Using (|4.10| ) we get from ( 4.11|) 



Xi-3(A)xi-3(Ae-"*) = 1 + e™ A (4.12) 

where a = argA/| argA| and < | argA| < n. 

The formula ( |2.6[) can now be used directly to define the Laplace transform Xi^3( s ) °f 
the Joos function xi->3(A) with the integration contour C13 running around the negative half 



of the real axis of the A-plane. In this way xi-+3( s ) is defined by (2.6) as the holomorphic 
function in the half-plane 3fts < 0. 




Fig. 12. The "first sheets" singularities of <b\ (£, s) for the harmonic potential 



The analytic structure of xi_>3(s) can, however, be best handled if we consider log* xi_»3(s) 
rather than the function itself (T^]. Namely, we have: 



log* xW*) = ^~. I e~ 2Xs log xi^WdX 



(4.13) 



and we can use ( 4,12| ) to calculate ( 4.13 ) exactly. (Note, that xi^3(A) does not vanish in the 
A-plane cut along the negative half of the real axis). Using ( [4.12D we have: 
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log*xi^3(s) = ^ J e- 2As log(l + e™ A )dA (4.14) 
C«(A ) 

+ i / e- 2As log(l + e-^ A )dA + ^ | e- 2As log X i^3(-A)dA 

C d (X ) C 13 (A ) 

where Ci3(A) is one of the contours C13 crossing the real axis at Ao > and C"(Ao), C d (Xo) 
are parts of it lying above and below of the real axis correspondingly. 




Fig. 13. The "first sheets" singularities of <b{ (£, s) for the harmonic potential 



Performing the integrations in the first two integrals in ( 4. (by expanding the loga- 
rithms) and changing A into -A in the third one we get: 



log*Xl-*300 



(-1) 



n+1 



n<l 



11 



_ inn 
2 



-2mn\o 



S + 



-2\ s 



(4.15) 



1 

'27d 



e 2Xs log xi- 



C"(Ao) 



i(X)d\ 



where C"(Ao) is the contour encircling (anticlockwise) the point A = and starting and ending 
at the point A = — Ao of the real axis. Since the right hand side of ( |4.15 ) is independent of Ao 
it can be calculated at Ao — * 0. It can be shown (see Appendix 5) that the integral in ( [4.15| ) 
vanishes in this limit and therefore we finally get: 



log* X1-3OO 



(-1) 



n+l 



n<l 



' ~l a— 



1 

2i~s 



1 

2is 



1 



sin (2is) 



(4.16) 
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The result §~W{) was established essentially by Voros ]l9| but here it is obtained directly 
by the definition ( |4. 13 ) of the Laplace transform for logxi^3(A). The inverse transformation 
can be also performed to give the known expression for %i^3(A) pTsfl . 




Summarizing the above analyses one can see that despite the clear way of obtaining 
corresponding singularity patterns and the underlying structures of the Riemann surfaces 
both they become still more and more complicated with increasing q. The following main 
observations follow, however, from this analyses: 

1. The set S q+ \ of singularities corresponding to contains the set S q of these 
corresponding to $^ (£, s). 

2. The new singularities which belong to S q+ i \ S q are generated on the sheets originated 
by the singularities of S q ; the latter is true both on the £- and on the s-planes. 
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The following comment concerning the positions of the singularities themselves and their 
relation to the Feynman path integral is in order. From the above discussion it is seen that 
these positions are determined by the values of the classical action the latter takes on along 
suitable classical paths corresponding to the case considered. The paths are real as well 
as complex (i.e. they are real or complex solutions to the classical equation of motion). 
They contribute to calculated quantities $j (£, s), q > 0, in a hierarchical way described 
above so that a path with greater absolute value of the real part of the corresponding action 
contributes to the later term (£, s) of the topological expansion. In this way the latter 
expansion reflects its close relation to the semiclassical expansion based on the Feynman path 
integral and the saddle-point technique as well as it confirms the role of complex classical 
paths in such calculations [15, [l6|, 17, [T§|1. 



5 An application: the connection problem 

The connection problem is an old problem of the JWKB theory which in the context of 
the Balian-Bloch representation was considered first by Voros |l9| . We shall discuss again 
this problem within the framework of our formalism to show the equivalence of the solution 
it provides with the corresponding method used in our earlier papers (see, for example, 



[12], [22],|2f]). 

The main question is how the JWKB formula, being a good approximation to a given 
solution in some domain of the x-plane, should be changed (in order to remain still the good 
approximation of the solution) when the solution is continued analytically to another domain 
of the x-plane. The problem can be solved in many different ways depending on the type of 



the considered solutions (see, for example, pO, 21, E3, pq[). In particular, it can be solved 
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with the aid of fundamental solutions (see |l2[ for the relevant procedure in an application 
to matrix element evaluations in JWKB approximation). 




Fig. 16. The Borel plane singularities corresponding to Xi (£) s ) 

In the framework of the Balian-Bloch representation the solution of the problem is the 
following. Consider the fundamental solution to ( p.l[) given by (^)- (|2.4|) and continued 
along a path 71" to sector 2 (see Fig. 2). As it follows from the previous section analysis 
continuing analytically along the considered path we can not meet singularities above the 
corresponding path 7i(£) in the £-plane and, therefore, the singularity pattern of Xi(£,s) in 
the s-plane looks like in Fig. 16a i.e. there are only two cuts on the relevant sheet. Since 
the integration along C is limited only to lie in the left half-plane we can deform it freely 
in this half-plane to the position C3, for example (see Fig. 16c). But doing this we have to 
integrate also along the cuts starting at the points s = £ — Ci and s = £ — £ii respectively. 
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Thus, A) is represented in this way as the following sum: 

A) = 2cri(£)e- A « / e 2Xs xi(Z,s)ds (5.1) 
•/ c 

2q~\{i)e-^(f +1+1 )e 2Xs xi(^s)ds 
\JCi JC 2 J Co, J 

where £ = q^dy. It should be noticed now that each term of the sum in (5.1) is a solution 
to the Schrodinger equation ( |2.1| ) (see, for example, ||]). It is also not difficult to see that 
each of the solutions generated by the integrations along C\ and C2 is proportional to the 
fundamental solution defined in the sector 2 of Fig. 2, whilst the remaining third solutions 
generated by the integration along C3 - to ip3(£,, A) - the fundamental solution defined in 
the sector 3. (An easy way to establish these facts is to investigate the behaviour of these 
solutions when £ goes to 002 and 003 correspondingly (oo^ being the infinity point in the 
sector k)). In this way a linear combination of the fundamental solutions ip2 and 1P3 to form 
the solution i/>i(£, A) is realized simply by moving the contour C in the s-plane. 

The connection problem arises when V , i(£> A) is continued to the sector 3 by crossing the 
sector 2 i.e. along some non-canonical path 71" in Fig. 2. At the end of such a continuation 
the dominant character of the JWKB factor q~ l / A exp(— A£) is lost in favour of the amplitude 
factor xi(C) A) but the series ( |2.3j ) does not give then an easy answer to what actually happens 
when K£ — > 00 along such a path. (In fact Xl(£> A) behaves then as exp(2A£)). 

In the s-plane the analytic continuation just described results in a deformation of the 
contour C to the form shown in Fig. 16d (broken line). It follows obviously from the figure 
that the dominant contribution to A) comes now from the integration along C2 i.e. from 
the solution proportional to ^/>2(£>A) and this is the way by which the connection problem 
is solved within the framework of the Balian-Bloch representation. (Note, that both the 
solutions defined by the integrations along C\ and C3 are subdominant when A — > 00 and £ 
stays in the sector 3 or when — > 003 and A is fixed). 

It is easy to see further, that the linear combination in the RHS of ( |5.1| ) can be explicitly 
reconstructed with the aid of the canonical coefficients a^j^ (a^/j^ = lim^^oo^, , see 
@) as follows: 

ipi(£, A) = ay 2 ^ p tp2^, A) + ai/j,->2^p(£ ; A) (5.2) 
= "1/2-^2 (£, A) + ai/p_ >2 a p /2^3V'2(£, A) + ai/ p ^ 2 a p /3^2V'3(£, A) 

where the sequence of terms of the last sum in ( |5.2[ ) corresponds strictly to the sequence of 
integrations along Ci, C% and C3 in Q5.1| ). The first linear combination appears in ( |5.2| ) when 
the contour C\ is deformed to the position C p shown in Fig. 16b. 

It is also worthwhile to note that the formula ( ^.2[ ) giving us the continuation of - 0i(£, A) 
to the sector 3 along the noncanonical path 71" can be also used to obtain in a simple way 
the improved connection formula of Silverstone p2[ (see also the recent work of Froman and 
Froman |2^j ) with ^3 playing the role of the subdominant contribution. It is enough to this 
end to substitute each term in the sums in ( |5.2| ) by its corresponding JWKB approximation 
(i.e. none cumbersome Borel resummation used by Silverstone is necessary). 
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6 Exponential asymptotics 



The problem of the semiclassical expansions for physical quantities is strictly related to 
the problem of so called exponentially small contributions absent (by definition) when only 
the bare semiclassical expansions of these quantities are considered [§3], |25|, ^6], . The 
exponentially small contributions become important if the accuracy of the best semiclassical 
approximation is considered to be insufficient. There are however two aspects of this problem 
a difference between which was, according to our knowledge, not discussed properly. 

The first one appears when the corresponding semiclassical series if Borel summed does not 
reproduce correctly the quantity considered. A discrepancy has to be of course exponentially 
small not contributing to the semiclassical limit. It can however be incorporated into the 
resurgent quantity by choosing another integration path in the Borel 'plane' i.e. the original 
path has had to be chosen incorrectly if the quantity considered was to be Borel summable 
(see || for the corresponding discussion). If it is the case then the Borel summation along 
the correct path has to incorporate all the exponentially small contributions to the quantity 
considered. 

The second aspect appears when the Borel transform reproduces completely the quantity 
considered and the semiclassical series is used as a source of the best approximation. As it is 
well known (see for example [p9| ) the latter is obtained in this case by abbreviating the series 
on its least term (since the series is divergent) the order of which is proportional to actual 
value of A(7i _1 ) (in fact n should be equal to the integer part of A|so| where so is a singularity 
of the Borel function closest to the origin). The remainder (i.e. the difference between the 
quantity and its approximation) is then exponentially small quantity. 

To improve this approximation using still the semiclassical tools we have to be able to 
identify the exponential factor of the remainder and to multiply the last factor again by some 
optimal abbreviation of a new semiclassical expansion of the remainder. Next, we should be 
able to repeat this procedure to the remainder of the remainder constructing in this way still 
more and more accurate semiclassical approximation which includes as many exponentially 
small contributions as we need to make the approximation as good as we wish. In Appendix 
4 we show how to do it. According to the beginning of this discussion the exponentially 
small contributions obtained in this way have to be determind by the singularity structure of 
the corresponding Borel functions provided for example by the topological expansions. The 
results of Appendix 4 confirm these expectations. 

7 Exponential asymptotics of energy levels 

Using the approximation scheme which follows obviously from the topological expansion 
and from the results of App. 4 we shall determine in this section the way of obtaining the 
semiclassical exponential asymptotic for energy levels of the anharmonic osciallator corre- 
sponding to the potential V(x) = x 2 +x 4 with the Stokes graph shown in Fig. 17 and drawn 
for E > 0. Taking into account the symmetry of the potential we can write the quantization 
condition for the energy levels in the form 

exp Q £ yjv(x) - E±(\)dx ± i|) = X1 ^ 3 (^ ± (A), A) (7.1) 
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where the contour K is shown in Fig. 17 and "db" in (7.1) correspond to the even and odd 
parities of the levels respectively. 

Imx 




Fig. 17. The Stokes graph for the harmonic potential V(x) 



+ x A for E > 



Making the complex conjugation of both the sides of ( |7. l|) we get an alternative condition 
for the energy level quantization. Both the versions are important since they determine the 
normal sector of E(X) to be < |argA| < 3ir/2 for A sufficiently large ||. Because of 
this the semiclassical series of E(X) is, as we have shown in our earlier paper [|j (see also 
Hi [| 0' H' Borel summable to E(X) itself and the singularity structure of E(s) on its 
Borel plane is also determined by ( |7. 1|) and its complex conjugation. As it follows from Q7. 1| ) 
this structure is symmetric with respect to the real axis of the s-plane and E(X) can be 
recovered by integrating E(s) along the negative halfaxis. Of course to apply to the last 
integral the procedure of App. 4 we have to know a detailed distribution of singularities 
of E(s) on its Borel plane. But instead of this we can use (7^1) directly to establish the 
respective exponentially small contributions to E(X). Namely ordering these contributions 
according to their exponential smallness we can treat each such a contribution as a correction 
to its predecessors and use the Taylor series expansion to take into account the corresponding 
contribution. So we can write 



E(X) = E (X) + E^X) . . . 



(7.2) 



with £i (A) being polynomially dependent on A 1 and with the contributions £i(A), E2 (A),..., 
being each exponentially small with respect to their predecessors. Of course, Eq(X) is con- 



structed in the standard way (see |11] and App. 4) and, for a given A, it has some well defined 
numerical value. 

The corresponding Taylor expansion of xi^3(E(X), A) with respect to the dependence of 
the latter on the energy is the following 



X i^ 3 (£(A),A) 



Xi^ 3 (E (X), X) + 9xi^(Eo(X),X) (Ei{x) + + _ _ _ } + _ _ _ {7 3) 
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Since the quantities in Q7.2| ) are real we shall write the corresponding quantization con- 
dition in its real form, too, to get 



sin (±fy/v(x)-E±{\)dx - ^(iSf (A) + £±(A) + . . ■)/ ; = + • • • ) (7-4) 

V K K 



dx 

\jv{x) - Eq{\) 



^ ( X i^(E ± (X), A) + (Et(X) + E±(A) + ...) 



d X i^(E±(\),\) 



+ ... 



It is now clear that we can apply to the coefficient xi-*z{Eq(X), A) and its derivatives the 
procedure of App. 4 considering E$(X) as having well defined value so that the singularity 
structure of the corresponding Borel function xx-^3(E, s) is determined just by the value of 
E equal to Eq(X). 

Assuming the exponential asymptotics for xi->3(Eo(X), A) to be ordered in a way analo- 
gous to ([T.2j) we get for the first two terms of (|7.2|) 



A 



sin ( f. j K sJV{x) - E*(X)dx j f ft (X& 3 (^< A)- A) ) ( 7.5) 



and 



Ef(X) 



±^( X £ 3 (^o ± (A),A) 



dx 



li 



cos U§ \Jv(x) - E^(X)dxj T « ( 8X ^ ( // (A) ' A) 



(7.6) 



where Xi^ 3 (.Eo(A), A) (we shall suppress the parity indices as unimportant for our further 
considerations) is given by the respective number of the first terms of the series ( [2.5| ) (i.e. 
abreviated at its corresponding least term; note also that the integrations in fl2.5| ) go from 
ooi to 003 along the canonical path) whilst the exponential contribution x^ 3 (.Eo(A), A) is 
determined according to App. 4 by the singularities of Xi^(Eq(X),s) in its Borel plane 
closest to the origin. 

Applying now the approximations following from the topological expansion flAl.ll ) for 
Xx— >3(-£'o(A), s) we can write (keeping only the first two terms of this expansion) 



where 



Xi_ 3 (£q(A), a) = <&£ 3 (#o(A), a) + $\% 3 (E (X), s 



(2) 



(7.7) 



$?l 3 (E (X),s) = I 



003 

As J tf(f,Eb(A))de 
\ 001 



003 

^l 3 (E (X,s))= J dq J d^ + V ,Eo(X))oj^,E (X)) 



001 



(7.8) 
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oo 3 / £ 6+?A 

z = 4( S -t?) | Co((,E (\))d( + 8(s- V )[J - J \l>^.EM\))<K 

OOl \OOl OOl 

Assuming the same order of approximation for Xi( x i -^o(A), A) we can see that when x 
stays in the sector 3 as it is shown in Fig. 18a then its Borel plane looks as in Fig. 18b on 
which C\ is the path of the Borel integration to recover Xi( x > _Eb(A), A). The distribution of 
the singularities on the figure follows now from ( [7.7]) . Fig. 18c shows the Borel plane for 

C(E (\)) 

Xi^3(Eq(\), s) i.e. when x — ► 003. The singular points are (c = J \/V( x ) ~ E-o(\)dx, 

B(EoW) 

C(E (X)) 

~Cc, Cc - (a = J \/V(x) - E (X)dx and ( A - Co- Therefore, xi^3(E (X), A) can be 
A(Eo(X)) 

given as 



Xi-*(Eo(\),\) = 2Ay e 2Xs xi^(E (X),s)ds 
a 



(7.9) 



where the integration path C is shown in Fig. 18c. 
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Fig. 18. The singularity patterns of 
Xi(x, Eq(\), A) for the anharmonic oscilla- 
tor (figure (a)) and xi, (x, Eq(X), s) (figure 
(b)). Figure (c) shows the latter pattern 
for x(£) — > oo 3 



The path C differs from the one considered in App. 4 but this does not prevent us to 
apply the procedure of this appendix. Therefore according to the approximation (7/7) we 
have 

Xi^ 3 (E (X),X) = 2XJ e 2Xs ($?l 3 (E (X), a) + ^l 3 (E (X), sj) ds (7.10) 

c 

so that 

X fl 3 (E (X),X) = E ^F^I ($?l 3 (E (X),s) + ^i 3 (E (X),s)) | s=0 (7.11) 

where n = [| ACc|] • 

Using the formulae ( |A4.5[ ) and ([A4.6Q of App. 4 for Xi->3(Eq(X) , X) we get 



(0) (E a\ A , _ V (no + 1)! ^ (-ir^ m) (^o(A),0) 

j=C,-C,A-C,C-A ( ZA ) ^ m=0 ^ 

where C-C = ~Cc, Ca-c = Ca~ Cc, Cc-A = Cc ~ (a, m = [\X(a\] (with \( A \ determining the 
common distance of singularities of Kj(Eo(X),s) closest to the origin) and Kj's are given by 

1 f ^l 3 (E (X),Cj+t) 



kJE (X),s) = / dt ^7 "V, 3 1 (7.13) 



(t + + ^ln(l + ^)- S ) (t + Ci + ai + ^ln(l + ^) 

tours ifj surround the cuts originated by 
Cj's each shifted to the origin s = of the Borel plane. 



~ (2) 

where the contours Kj surround the cuts originated by the singularities of ^^(E^A), s) at 



28 



8 Summary 

In this paper we have found the representation for the Borel functions of the quantities 
relevant for the ID quantum mechanics. The representation takes the form of the topological 
expansion. This expansion provides us with an algorithm determining in a systematic way the 
singularity structure of the Borel plane for the relevant quantities and orders the appearing 
of the Borel plane singularity structures in a hierarchical way allowing for the formulation 
of the approximation scheme of the semiclassical calculations alternative to the other ones 




We have also remarked limitations of our method in the description of the proper nature 
of singularities of the quantities represented by the expansion. 

We have formulated also the scheme of the semiclassical approximations including the 
exponentially small contributions to the desired order of accuracy. It makes use of the Borel 
plane singularity structure in the most natural and effective way, particularly, if it is accom- 
panied by the topological expansion method of approximations of the Borel functions. 

We have demonstrated the action of both the methods considering some simple (but not 
quite trivial) examples of their applications in Sections 4-5 and 7. However, it was not our 
aim in this paper to perform some numerical tests of the method presented. Rather we have 
limited ourselves to test both the methods as theoretical tools for better understanding of the 
mutual relations between the semiclassical expansions, Borel plane singularity structure and 
the exponential asymptotics. For the latter goal both the expansions (i.e. the topological 
and the exponential ones) appeared to be very useful. Nevertheless, their test as a practical 
method of extended semiclassical approximations is certainly desired. 
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Appendix 1 



Al.l. The Laplace transforms Y n ;ri,...,r q {£>s) 
We shall determine below the Laplace transforms Y n ;ri,...,r q {^ s )> n > q > 0, as defined 
by ( |3.5| ). To begin with consider first the case r% = and q = 1. We have: 



(Al.l) 



71 (0 



7l(£n-l) 



r = l,...,n, n = l,2,... etc. 

The multiple integral in ( [Al.l ) can be rewritten further as follows: 
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; V(£,A)= J d^ r e 2X ^^ r )n r ^^ r )Y^ r ) (A1.2) 

71 (0 

r = l,...,n, n = 1,2, ... etc. 
where Y^} r .{^ r ) = Y n _ r; (£ r ) is defined by ( |3.3| ) and 

fir-iCCCr) = ((r-l)!)- 1 TO-n($ r )) r - 1 (A1.3) 
Making further in ( Al.2| ) a change £ r — > £ r — s of the integration variable we get finally: 



4V & A ) = J dse 2Xs Y$ (A1.4) 
c 

r = l,...,n, n = l,2, ... etc. 

where 

4V(^) = - s)nr_i(e,e- 8)Y£%.(Z - s) (A1.5) 

and the contour C runs from the infinity = — oo to the origin s = 0. Note, that because 
of > (by assumption) the contour C is independent of r and n and also, as it follows 
from (|A1.5| ), injr (£>s) is holomorphic for 5Rs < 9?£. 



Resoning in the completely similar way we get for Yn;r ir2 \£, s): 



Wink**) = - J + s)w(e 1 )n ri -i($,ei + s)n r2 _ ri _i(^i + s ,a)i; ( ° ) r2; (6iAi.6) 



The remaining Laplace transforms in,T'!i 2 , + i (£, s) and i^jn... r 29 (^j s), q = 1,2,... etc. 
can be defined recurrently as follows: 

Y™r 2 r 3 ...r 2q &s) = - J df] J ^fe + «+#r 1 -l((,6+«) (A1.7) 

£7(s) 7(£-«) 

^r 2 -n-l(?l + S, £l + »7)4-^-r3...r2,-ra (& + ??,??), 9 = 2, 3, ... 

and 

4™...r 2q+1 (e,^) = - / d^(e- fl +7 ? )i2 P1 -i(e,ei-a+'7) (ai.s) 



£7(5) 



^ / n-ri;r 2 -ri...7-2„+i-r :L (£i s + ??) 7 ?) q— 1,2,... 
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Every of them is holomorphic in £ and s for > and 5Rs < The contour C(s) in 
(Al.7)-(A1.8) starts at the point s with "Sis < and ends at s = 0. 



A 1.2. Topological expansion 

A further step we can do is to fix q and to take sums with respect to n, ri, ...,r2 q . It 
can be done as follows. First, we consider rather A)/(2A) than A) itself. Next 

we note that to each term — (— 2\)~ n ~ 1 Yn?r 1 ...r q {£,, s) there corresponds the following Laplace 
transform: 



where the star means the convolution of the factors. 
The sums we are now looking for are the following: 



= E ( 1} ' J' ' > *f,&...r«(6«) (ALIO) 

so that the series: 



ri+r 2 +...+r q 



Til 

l<ri<...<r q <n 



$i(H,s) = J2^i\ts) (ALU) 

(its convergence is discussed below) represents a function s) such that cM>i(£, s)/ds is 
the Laplace transform of Xl(£> A). 

The sums in ( |A1.10| ) can be performed explicitly to give: 



^i 29) (C,s)= y dr/i | d m ... J drj q J dt x jd&... J di q 



C(s) C{ m ) C( Vq -i) 



oo oo oo 



^(6 + +%)••••• + 77 g )a)(^)(2s - 2?] 



U 3 

Z 2<? 



z 2g = 4(s - mMO + 8(s - m) £ (fi(£ P + ^+1) - + »&>)) . 

p=l 



7fo+i = 0, g = l,2,... (A1.12) 



^»i 2<?+1) (e, ^) = / <% • • • / dr, q+1 u>(£ -m + m) J #!•••/ *9 



C( S ) C( % ) 
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z 2q +i = 4(s - »7i)fi(0 + 8(s - 771) £ (fl(fp + r/p+a) - 0(£ p + 7^+1)) 

p=0 

£0 = £, = 0, g = 0,1,2, ... 



The functions I q (x), q > 0, in (A1.12) are the modified Bessel functions (of the first kind, 
see [31] p. 5, formula (12)). The results ( A1.12 ) have been obtained from ( |A1.10| ) by using 



repeatedly the following sum rule [31 



(A1.13) 



valid for any v. 

The formulae ( A1.12| ) provide us with the general forms of $1 (£, s)'s. However, for the 
singularity analysis of the latter the more convenient representation for them is the following 
recurrent one: 



$>f q+2 \i,s) = - J drj j drj J dr }1 u^ 1 )u^ l -r,'){2s-2ri) 

C(s) C(n) T/(ri) 



^ h (V4( a - y)(n(Q - 2n(gi) + n(& - if 
y /4(s-ri)(n(0-2n(£ 1 ) + n(z 1 -Tf)) 

$ { i q+1 \t,s) = - J drj J drfu>(e-rf) (ALU) 

C(s) d( v ) 



9 = 0,1,2, 



where ^\^s) is given by ( [Al.l2j ). 

Note that ( [A1.14 ) can be obtained from ( A1.12j ) and vice versa by applying the following 
relations: 



dxI m (^)I m (J (3(1 - x))(ax^ m ((3(l - x))h = 2 tt 7 Jm |^ | (A1.15) 



(s- V ) n 



1 



(k-l)\(n-k)\ 



(VaT?) m+n+1 
drf{ 8 -rf) kr - 1 {rf - nT~ k 



A1.3. Analytic properties of the functions (£, s) 
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Since each of the functions I q {z\^ 2 ) / ' zj/ 2 , q > 0, is an entire function of its argument 

then it follows from ( |A1.12| ) that possible singularities of $j (£,s) are generated by the 
(known) singularities of the functions 0^(77) and 0(77) and their integrations present in ( A1.12| ). 
However, it can be easily checked that the conditions: 



K£ > and < (A1.16) 



determine the domain where the integrands in ( [A1.12 ) are holomorphic. Therefore, this 
is also the domain of holomorphicity of (£,s) since all the integration paths in ( |A1.1| ), 



( A1.12 ) can be chosen to lie completely in this domain. 

Let us note, however, that as it follows from ( A1.12| ) each s), q > 1 can be con- 



tinued analytically from the domain ( |A1,16| ) to any point £ of the £-Riemann surface of 
$j (£, s ) (obtained for fixed s) if the distribution of branch points of cD(£) along a path of the 
corresponding analytical continuation is such that the distance of any of them from the path 
is greater than This statement is the direct conclusion from the corresponding formulae 
in ( [Al,12| ) since all the integrations on the £-Riemann surface present there are performed 
inside a strip no wider than \s\. 

Let us note further that for the polynomial potentials the branch points of £o(rj) are 
isolated and on each sheet of the £-Riemann surface of ib{rj) their numbers are finite. The 
distances between them on each sheet are nothing but the corresponding distances between 
turning points measured by the action. Therefore, there is the smallest distance d among 
them. If we take, therefore, s in ( |A1.12 ) such that \s\ < d' < d/2 then we can penetrate by 



paths of the analytical continuations the whole £-Riemann surface of <!>(£) if the former is 
deprived all the circular vicinities of radius d", d' < d" < d/2, centered at each branch point 
of the surface. We shall denote the corresponding part of the £-Riemann surface as R(d"). 
Consider now a question of convergence of the series in (Al.ll). We shall show below 



that the series is convergent absolutely and uniformly in the domain R(d"). It means that 



the series (Al.ll) determines ^(^s) as the holomorphic function in these domains. 

To this end let us note that if \s\ is chosen to satisfy the condition \s\ < d' < d" all the 
integration paths 7(£— f]i) can be deformed then to lie inside an infinite strip 5(£, s) bounded 
by the paths 7(£) and 7(£ — s) so having the width \s\ with the one end of the strip being 
placed at the infinity 00 1 and the other one being a segment (£, £ — s). The latter bound 
can be chosen as such because the path C(s) can be deformed to a segment (with its ends 
anchored at the origin and at s). Introducing now the following functions: 



H(£r,f7i) = hmsup |w(£ r + r?x) | (A1.17) 



\p\(t,Vl) = J \d£r\\ti\(€r,V) 



we have: 



71 (0 



m r + rj)\<p^, m ), rj G (7(771), £ r e 7i(0 (A1.18) 
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and for q large enough: 



\z q \ < 8(q + l)\s-Tn\p(Z,rn) 



(A1.19) 



so that: 



1 



|2% 2 i q \ < - 7 ex P (2|s- mlpfom)) 



\s\ 



2q 



1*1 (£> s )l < ({2qy.ql(q - I)])' 1 J dxx q -\\s\ -x) 



/ \* 

y I^P| 2 (^,m) I <'XT>(2(M -./•)/>(£. /; L )) 

\h(£-'7l) 



where x = \r]\\. 

Introducing yet: 



-i 



1*1 



< ((2 9 + l)!( g !) 2 J y d.v.rn\*\-.r)^\^\( lh <h) 




|dr/||u;| 2 (r/,r/i) exp(2(|s| - J71)) 



\7i(?-'?i) 



(A1.20) 



= hmsup|o}|(C,r?i) 

r?ieC(s) 

= limsup/5(f,r/i) 

«ii£C( s ) 



s) = limsup / 

we obtain finally for q — > 00: 



|rfry||cZr|(77, 771) 



(A1.21) 



$ i 29) ^ s )l<7frr7^ 9 (e^)e 



|3g 



2<7ft e \„2|s|p(£, S ) 



*? q+1) (Z,s)\< 



\3q+2 



(3q + 2)\ql 



(3q)\q\ 

Q 2 «(£, S )|u,|(e, S )e 2 H^> s ) 



(A1.22) 



The bounds ( |A1.22| ) show clearly that the series (ALII) is convergent in the assumed 
domain R(d") since Q(£,s), s) and | (77, s) are finite there. 
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Appendix 2 



If x p is a simple zero of q(x) then the point £ p = £(xo, x p ) is the branch point for the function 
u)(£) defined by ( f^ ) which can be expounded around the point £ p into the following series: 



s(o= E ^(e P )(^-^) 2fc/3 (A2.1) 

fc>-3 



The coefficients cjfc(£ p ) in (A2.1) are defined by the identity: u>(£(xq,x)) = tu(x)q 2 (#) and 
the following expansions of £(xo,x) and u;(a;)g~2 (xc) (see (|2.4j)) around x p : 



£(x ,x) -£ p = E^ fc ( x p)( 



3 



k>0 

and (A2.2) 



io(x)q 2(x) = ^2uj k (x p )(x - x„) k 3 



fc>0 

In particular, the coefficient at the most singular term in (A2.1) cu_3 = —5/36 i.e. it is 
potential independent. It depends, however, on the multiplicity of zero of q{x) at x p , namely, 
lj_3 = —n{n + 4)/[4(n + 2) 2 ] for the n-fold zero. 

Appendix 3 

We establish here the Riemann surface structure of $j (£, s) for the linear and harmonic 
potentials. 

3. 1 The linear potential 

According to Section 4.1 the Riemann surface structure of $^ (£, s) for this case is de- 
termined by 



*s"«-)=-l/^f^ \ 7 (A3 ' 1) 

C(s) 4 £ 



From ( |A3.1 ) it follows that its subintegral function is singlar at £ = rj and at £ = where 



- 1 7 - 1 1 

it behaves as e ± ^~ r '- ) 3 (£— 77)" and e 1 *^ 3 £2 respectively. The ^-integration generates only a 
singularity at s = £ (by the EP-mechanism) leaving the singularity at £ = nad its character 
unchanged. Therefore, assuming £ to be continued to the sector 2 (along the canonical path) 
the corresponding first sheets of the Riemann surface look as in Fig. 10a, b. 

~ (2) 

Consider now &\ (£, s). Its Riemann surface structure is defined by: 
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C(s) 7(0 

5, ,/l 2 2 



As it follows from (A3. 2) the subintegral function is singular at £i = rj, £ = and at 

£i = behaving there as e^ 1- ^ 5 (£i — 1 £2 and e^ 1 ^ 1 respectively. 

The £i-integration in ( |A3.2|) generates the EP-singularities at £ = 77 and at £ = but also 
the P-singularity at 77 = (when the singularity at £1 = 77 move around the end point of 7(£) 
clockwise pinching the latter against the singular point £1 = 0). 

The final 77-integration in ( |A3.2j ) generates the EP-singularity at s = £ and s = and the 
P-singularity at £ = 0. Therefore, the 'closest' singularities of (£)-s) are the following: 

£ = 0, £ = 77, s = (A3.3) 

Let us make a general note that the EP-mechanism repeats the distribution of the branch 
points and cuts whilst the P-one generates new branch points on the Riemann surfaces 
obtained by the EP-mechanism always however enforcing specific ways of moving around the 
singularities generated by the EP-mechanism. 

All the singularities in ( A3.3| ) are the root branch points (of the forth order) accompanied 



by essential singularities as we have mentioned above. The singularity at s = 0, however, 
to be reached needs to round the branch point at s = £ moving clockwise i.e. it lies on the 
sheet opened by the latter branch point. This singularity is therefore a consequence of the 
singularity at £ = 0. 

A similar note concerns the singularity at £ = 0. In fact there are two such singularities 
the one on the sheet shown in Fig. 10a (arising by the EP-mechanism) and the second one at 
the sheet opened by the branch point at £ = 77 i.e. to reach it one needs to round this point 
clockwise. 

One can conclude therefore that the P-mechanism applied once has generated a singularity 
at s = and applied twice has generated a new singularity at £ = (on a different sheet) 
from the old one. It is clear that this mechanism will proliferate the last singularity on all 
the sheets of the Riemann surface of $i(£, s) except the sheet we have started with on which 
the point s = is regular for $!(£, s). 

Now we can use the formulae flAl.14 ) to prove the form of the first sheet of the Riemann 



surface as shown on Fig. 10a, b. Namely, assuming for &^ q \l;, s) the form of this sheet shown 
in the last figure we deduce that it remains unchanged for $^ 2<7+2 ^(£, s) whilst it is deprived 
of the singularity at s = for &^ q+1 \l;, s). 

Indeed, consider the subintegral function in the first of the formulae ( Al.l4j ) defining 
3>^ 2<?+2 ^(£, s). It has singularities at the following points: 



£ = 0, 6 = 0, £i-t/ = 0, £1-77 = 0, r,- V ' = (A3.4) 
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shown on Fig. 19a for the £i-Riemann surface (when the rest of the variables are fixed). 



Im|, 



Im r| 



Reti' 



Fig. 19. The "first sheets" singularity pattern of the subintegral function defining 



(£, s) for the linear potential before the ^-integration (a) and after it (b) 



The ^-integration provides us with the EP-singularities at £ = 0, £ — 77 = and £ — n' = 
and with the P-ones at 77 = 0, rj' = (when the point £ is rounded by 77 and 7/ clockwise to 
touch the point £ = by the latter) and at 77 = rj 1 (when the point 77(7/) rounds £ clockwise 
(anticlockwise) to touch 7/' (7/)). We get in this way the singularity pattern before the 77'- 
integration shown in Fig. 19b where the singularities at 7/ = and at 7/ = 7/ are screened by 
the £-cut and to get them one has to go around £ clockwise or anticlockwise respectively. 

The T/'-integration therefore does not do much now providing us with the singularity at 
7/ = £ and at 77 = (by the EP-mechanism) and at £ = (by the P-one) with the latter 
singularity placed on a sheet opened by the singularity at £ = 77. 

The final ^-integration repeat only the singularity pattern described just above so we are 
left with the distribution of the singularities as shown in Fig. 10a, b. 

Consider now the second formula (A1.14). There is no the ^-integration and therefore 
the singularity at rj' = is not generated and the other singularities at £ = 0, 77 = and s = 
can not be generated as well by the further 77'- and 77-integrations. Besides the generation 
of the singularities at £ = s goes exactly in the same way so that the final picture of the 
corresponding Rieman surface is the same as in Fig. 10a, b except of missing of the respective 
singularities at £ = and s = on the lower sheets. 



3.2 The harmonic potential 



We assume here 7i(£) to be continued canonically to the sector 3 of Fig. 11 and we put 
= I-i Vx 2 + ldx = £. Nor a)(£) nor S7(£)are now simple functions of £. a)(£) is periodic 
(with its period 2£ acting between different sheets of the infinitely sheeted Riemann surface 
on which this function is defined) whilst is not. 



37 





Im r| 



Rei; 



Re n 



C(n) 



Fig. 20. The singularities of the subintegral function in (|4.1| ) defining $1 (£, s) 
for the harmonic potential 



Consider however again s) as given by (4.1). The closest singularities of the subin- 

tegral function are shown in Fig. 20a, b i.e. they are: 



C = 0, £-C = 0, £-77 = 0, e-77-C = 



(A3.5) 



Therefore the ^-integration in (4.1) provides us with the following singularities of (£, s) 
shown in Fig. 12a, b: 



e = o, e - c = o, s-8 = o, £-c- s = o 



(A3.6) 



i.e. none the P-singularity is generated. 

Consider next 3^ (£, s). According to (f4.3|) singularities of the subintegral function in 
this formula are now the following: 



e = o, e - c = o, 6 = o, d-c = o 

6. - ?7, £i - C - V = 



(A3.7) 



and we can use Fig. 20 to show the corresponding situation situation for the £i- and 77- 
dependance by making the substitution £ — > £1 in the figure. 
The ^-integration generates the following singularities: 



C = 0, e-C = 0, £-r7 = 0, £-£-^ = 



(A3.* 
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by the EP-mechanism and 



77 = 0, 77 — c = 0, ^ + ( = 



(A3.9) 



by the P-mechanism. The latter singularities lie on the 'lower' sheets of the r/-Riemann 
surface. 

Finally, integrating in ( |4.3[) over 77 we generate singularities of s) at 



£_ s = 0, f-C-s = 0, s-C = 0, s + c = o 
by the EP-mechanism and at 



(A3.10) 



e + c = 0, e = o, e-c = o, e-2c = o 



(A3.ll) 



by the P-mechanism. The proper distribution of these singularities is shown on Fig. 13. 
Now we can proceed inductively assuming for $1 s) the singularity pattern shown 



in Fig. 14a, b where is the integration path in the formulae ( |A1.12| )-(A1.14) and C the 
corresponding path to recover from (£, A) by the Borel transformation (at s = (£, s) 
is then regular). 

Taking into account the second of the formulae ( |A1.14| ) we see that the singularities of 
the subintegral function are determined mostly by its factor ^^ q \^ — rf, rj — rf) according to 
which and Fig. 21a these singularity are at the points 



£-77 = 0, £-r,-C = 0, £-V-K = 0, k = -(2q-l),...,2q 
rf - rf + k( = 0, k = —(2q - 1), (2q - 1), k^O 
shown for the case of the corresponding r/'-Riemann surface on Fig. 21b. 



(A3.12) 



b) 




Re I, 




.ccn)_ 



Im if 



Re 11' 
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5+(2q-l)S 



Im n 



^-2 



Re r| 



Fig. 21. The singularity structure of 

^l. 2 ^(£ — v')V~ v') f° r the harmonic poten- 
tial in the £-plane (a), in the //-plane (b) 
and in the //-plane after the //-integration 

(c) 



Therefore, making the r/'-integration we obtain the '//-plane' singularity pattern shown 
in Fig. 21c on which the ^-dependent singularities are created by the EP-mechanism whilst 
the two fixed ones on the imaginary axis by the P-mechanism. Other singularities generated 
in the last way appear on the lower sheets originated by the two singularities at r/ = £ and 

v = C- 

The successive //-integration does change nothing in the s-variable singularity pattern (in 
comparison with this on Fig. 21c) so providing us finally with its form shown in Fig. 15b 
but seriously changes the original pattern of Fig. 14a. Namely, the EP-mechanism generates 
the ^-singularities at the points £ = s — (2q — 1)£, s — (2q — 2)£, s — £, s, s + £, s + 2qQ 
and by the P-mechanism at the points £ = — 2g£, 0, (2q + 1)£. As the final result 

we have for ^^ q+1 \^, s) the picture of Fig. 15 for its both types of singularities. 

The corresponding analysis of the case ${ 2g+2) (£, s) is s till a little bit more tedious but 
nevertheless direct due to the first of the formulae ( |A1.14| ). The valid singularities of the 
subintegral function in this case are 



£1 = 0, £1 - 77 = 0, 

£l-7, = 0, £1-77- £ = 0, £1 - rf - K = 0, k = — (2<7 - 1), 2q (A3.13) 

rf - 7/ + k( = 0, k = —(2q - 1), (2q - 1), k^0 

The first £i-integration is performed on the sheet shown in Fig. 22a. By the EP- and 
P-mechanism it generates 77- and //-singularities. Limiting ourselves to collect only these 
singularities which appear on the //-sheet on which the result of this £i-integration is regular 
at 77' = we arrive at the pattern shown in Fig. 22b. 
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Fig. 22. The subintegral singulari- 

ties in ( |A1.14D defining $f 9+2) (£, s) in the 
£i-plane (a), in the r?-plane after the £1- 
integration in (A1.14) (b) and in the 77- 
plane (c) after the ^'-integration 



The successive r/-integration leads us to the '77-plane' pattern shown in Fig. 22c where 
the two ^-dependent singularities were produced by the EP-mechanism whilst the fixed ones 
by the EP- and P-mechanisms simultaneously with the exception of the highest two ones at 
77 = — (2q + 1)C and 77 = (2q + 1)£ which are generated by the P-mechanism only. 

The final 77- integration in ( |A1.14 ) provide us with a pattern analogous to the one of Fig. 
14a by the EP- and P-mechanism and with the pattern of Fig. 14b by the EP-mechanism 
when q is substituted by q + 1 . 



Appendix 4 

We describe here a procedure allowing us to construct in a systematic way the opti- 
mum semiclassical representation for the Borel summable quantity including both the main 
contribution coming from the semiclassical series abbreviated at its least term and the cor- 
responding exponential contributions of an arbitrary order. In its finite form the procedure 
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provides us with the exact formula for the quantity considered. However, if continued in- 
finitely the procedure give rise to the question of convergence of the infinite functional series 
we get by it. 

To this goal we shall consider the basic quantity given by the formula ( |2.8| ). Integrating 
in it by parts we get 



/ 



Xi(£,A)=2A 



E 

k=0 



\n+l 



(2\y 



(2\y 



n+l 



,2As 



(A4.1) 



c 



According to the well known prescribtion (which can be easily justified by the analysis 
similar to the one performed below) we should put n = no = [A|so|] in (A4.1) where [x] means 
the integer part of x and so is a singularity of s) closest to the origin. Next we should 
extract from the integral the exponentially small factor and finally continue the procedure to 
the remaining Borel integral. This could be done in the following way. 

First the no + 1 th derivative of Xi(£, s ) can be given the form 



X { r +1 \t s) = ( - 1)n Qr +1)! J* frit, * + t)t-^dt (A4.2) 

with the integration contour K in ( A4.2| ) surrounding anticlockwise the negative half axis of 
the Borel plane (see Fig. 5). 

As it follows from Fig. 2 sq = £ — £i = £ (since £i = 0, see Fig. 4). Deforming the contour 
K to surround the cuts generated by the points £ — (i and £ — £2 of Fig. 5 (for £ chosen as 
in the figure they are the unique cuts visible in these positions) and shifting the integration 
variable in the corresponding intagrals we get 



xi" +1) (^) = ( l)no+ { l^ + l)l i: J K m-Q+m-Q-s+t)-^dt (A4.3) 



where the contours Kj surround (anticlockwise) the cuts whose origins are at the point t = 0. 

Further substituting ( |A4.3|) to (|A4.l| ) and changing both the order of integrations in 
( A4.3| ) and the integration variables themselves we get as a result of these calculations 



Xifo A) = E y^*i%0) - g {2X )n \l%n J e 2 ^,s)ds (AAA) 

h — J 1 C 

where 

1 f M Xl(.te-<3+t) 



K, f,5 = / dt ; t 

' 2m J (l + j= 



no 



(t+e-o + ^in(i + ^-)- s ) (t+e-o + T + i? ln (i + ^cj)-« 

3 = 1, 2 
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(A4.5) 



and where the contours Kj run again around the cuts anchored at t = 0. 

The form ( A4.5| ) for k's allows us to continue the procedure of getting the asymptotic 



series expansions for the integrals in ( A4.4| ) and next to abbreviate the series at their least 



terms. The latter are to be determind by the singularities generated by the i-integrals in the 
s-plane (as a result of the pinch mechanism) closest to the origin of the plane. It is easy to 
see that among possible candidates for the latter are the singularities at s = £, £ + no/A for 
s) and the ones at s = £ — C2, £ + no/A — C2 for K 2(£ 5 s) (all the singularities are generated 
by the P-mechanism at t = 0). However, the integrations in ( |A4.5 ) along the corresponding 



cuts open possibilities for new singularities to appear generated by the t-singularities shared 
by the cuts. These possibillities still enrich the variety of singularities which have to be taken 
into account in choosing the one closest to the origin of the s-plane. 

Therefore, to construct the representation (A.4.1) for each of the two integrals in ( A4.4[) 



we have to choose from the singularities corresponding to each k the ones which are closest 
to the origin. When these choices are done the procedure described above can be repeated. 
Let us call Kj, j = 1, 2, defined by ( A4.5| ) the first generation family considering kq(£, s) 



Xi(£, s) as the zeroth generation one. It is clear that the general form of the optimum 
semiclassical representation for xi(£> A) is the following 



m=0 



Xi(£,A)=2A^L (2A)-+i 



2aV(-d* y n + y ( - 1)m ^..,i fc (^) (A46) 



k=l ji,-,jkl= 



Jl>— ,3p+ 



ii t 1 (2A)^- 1 ^(£-0 



where jo = and Kj 1 ...j p+1 's constitute the p+ 1 generation family. The latter is constructed 
from the p th one (with (j p as its singular points and with nj p /X = |CvJ being the singularity 



closest to the origin) according to the formulae ( A4.2| )-( [A4.5| ). 



It is important to stress that ( |A4.6| ) is exact and its RHS becomes an approximation to 
the left one only when the last sum of the RHS containing the integrals is rejected. 

Appendix 5 



We shall show below that the last term on the right hand side sum in ( 4.1 5| ) has to vanish 



when Ao — > 0. To this end let us note that we can rewrite the integral present in this term in 
the following way: 



J exp(2As) log xl -» 3(A)d\ = J [exp(2As) log *1 -» 3(A) 



C'(Ao) Ci 

2 



exp(-2As) log xl -» 3(-A)]dA = J [exp(2As) - exp(-2As)] log xl -» 3(A)dA (A5.1) 



Ci 

- 
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+ J exp(-2As)[l + exp(-2vniA)]c2A + J exp(-2As)[l exp(2vra)]dA 

2 2 

where Ci is the half-circle of radius Ao lying in the right half of the A-plane, and Ci and Cf 

2 2 2 

are the corresponding upper and lower halves of Ci. We have also made use of the relations 
( 4.11 ) to obtain the final form of ( A3.1| ), It follows now from ( 4.10| ) that we have: 



lim xl -» 3(A) = \/2for \arg\\ < vr (A5.2) 

A — >0 

Therefore, we can conclude that both log\xi—>3(X)\ and argxi^z(X) are bounded in the 



halfplane 5RA > 0. The vanishing of all the integrals in ( A5.1 ) when Ao — > follows now 
directly from the last conclusion. 
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